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Ich möchte in diesem Bericht algorithmische Methoden vorstellen, die im wesentlichen in 
diesem Jahrzehnt Einzug in die Computeralgebra gefunden haben. Die hauptsächlichen Ideen 
gehen auf Stanley und Zeilberger ||3^-[^ zurück, vgl. die Beschreibung und haben 



ihre Wurzeln teilweise bereits im letzten Jahrhundert (siehe z. B. 0-0), gerieten aber auf 
Grund der Komplexität der auftretenden Algorithmen wieder in Vergessenheit. 
Eine der Kernfragen kann hierbei so formuliert zu werden: Worin liegt der wesentliche Unterschied 
zwischen der Exponentialfunktion f{x) = und beispielsweise der Funktion g[x) = + 
|x|/10^°°°, der dazu führt, daß alle Mathematiker die Exponentialfunktion / als elementare 
Funktion auffassen und nicht g, obwohl sich diese beiden Funktionen auf einem großen Bereich 
der reellen Achse numerisch kaum unterscheiden? 

Oder ein Beispiel aus der diskreten Mathematik: Warum erhält die Fakultätsfunktion a„ = n! 
den Vorzug gegenüber &„ = n! + oder irgendeiner anderen diskreten Funktion? 

Obwohl die Beispiele die bekanntesten stetigen bzw. diskreten transzendenten (nicht- 
algebraischen) Funktionen betreffen, ist die Antwort auf diese Fragen rein algebraischer 
Natur: Die Exponentialfunktion / ist nämlich charakterisiert durch jede der folgenden algebraischen 
Eigenschaften: 

1. / ist stetig, und für alle x, y gilt f{x + y) = f\x) ■ f{y); 

2. / ist differenzierbar, und f'{x) = f{x) sowie /(O) = 1; 

oo 

3. / G C°°, f{x) = J2 o-nX^ mit ao = 1, und für alle n > gilt {n + 1) ün+i = ün] 

n=0 

und die Fakultätsfunktion wird charakterisiert durch jede der folgenden algebraischen Eigenschaften: 

4. ÜQ = 1, und für alle n > gilt a„+i = {n + 1) a„; 

oo 

5. die erzeugende Funktion f{x) = J2 dnX^ erfüllt die Differentialgleichung x'^f'{x) + 

n=0 

{x — l)f{x) + 1 = mit der Anfangsbedingung /(O) = 1. 

Während mir keine Methode bekannt ist, Funktionalgleichungen wie Eigenschaft (1.) in der 
Computeralgebra zur Darstellung transzendenter Funktionen einzusetzen, will ich aufzeigen, 
in welcher Form die anderen Eigenschaften dazu geeignet sind. 

Man beachte, daß die erzeugende Funktion der Fakultät nur am Ursprung konvergiert, also als 
formale Reihe anzusehen ist. D. h. insbesondere, daß man eine geschlossene Darstellung der 
erzeugenden Funktion nicht angeben kann. (Diese ist zumindest nicht eindeutig: Jede Funktion 
der Funktionenschar — ^-^(Ei(— l/a;)+c) (s. [|I|, Kapitel 5) stellt die erzeugende Funktion der 
Fakultät für x > asymptotisch dar.) Darauf kommt es aber auch gar nicht an: Statt mit der 
erzeugenden Funktion selbst arbeitet man ohnehin viel besser mit ihrer Differentialgleichung, 
welche rein algebraisch ist! Genau dasselbe trifft auf die Exponentialfunktion und die Fakultät 



1 



selbst zu: Statt mit diesen transzendenten Objekten arbeite man mit den dazu äquivalenten 
algebraischen Differential- bzw. Rekursionsgleichungen (vgl. ^j). 

Die gegebenen Eigenschaften sind Strukturaussagen für die betreffenden Funktionen. Bei 
der geringsten Änderung geht diese Struktur verloren. Beispielsweise kann die Funktion g{x) = 
+ |a;|/10^°°'' durch keine den Eigenschaften (l.)-(3.) analoge Vorschrift charakterisiert 
werden. Dagegen ist die Funktion h{x) = e^+x/10^^^^ beispielsweise durch die Differentialgleichung 
(x — 1) h"{x)—xh'{x)+h{x) = mit den Anfangsbedingungen h{0) = 1 und h'{0) = 1 + 10~^°°° 
gegeben. 

Das Besondere (und Gemeinsame) an Exponential- bzw. Fakultätsfunktion besteht also darin, 

daß diese eine homogene lineare Differentialgleichung und jene eine homogene lineare Rekursionsgleichung 

erfüllt, wobei Differentialgleichung sowie Rekursionsgleichung Polynomkoeffizienten haben 

und beide erster Ordnung sind. 

Verallgemeinern wir diesen Sachverhalt nun zunächst auf stetige Funktionen einer Variablen 
und nennen eine Funktion f{x) holonom, falls sie eine homogene lineare Differentialgleichung 
mit Polynomkoeffizienten in x erfüllt. Stanley püf zeigte, daß Summe und Produkt holonomer 



Funktionen sowie die Komposition mit algebraischen Funktionen wieder holonome Funktionen 
liefern. Beke hat bereits vor 100 Jahren Algorithmen beschrieben, mit welchen die 

Differentialgleichung für Summe bzw. Produkt von / und g aus den Differentialgleichungen 
für / und g bestimmt werden können! 

Analog nennt man eine diskrete Funktion a„ holonom, falls sie eine homogene lineare Rekursionsgleichung 
mit Polynomkoeffizienten in n erfüllt. Summe und Produkt holonomer diskreter Funktionen 
sind wieder holonom, und es gibt Algorithmen zur Berechnung der entsprechenden Rekursionen 
(s. 
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Was haben wir hiermit nun gewonnen? Ignorieren wir einmal, daß e^, sinx, cosx, arctanx, 

aresin x etc. transzendente Funktionen sind, und stellen wir lediglich ihre holonomen Differentialgleichungen 

/' = /, /" = -/, /" = -/, (1 + + 2xf' = 0, {x^ - l)f" + xf = etc. in Rechnung, 

so können wir nun aus diesen Differentialgleichungen mit reiner Polynomarithmetik (wir 

brauchen eigentlich nur lineare Algebra, vgl. [^, ||2^) holonome Differentialgleichungen für 



Summen und Produkte solcher Funktionen, beispielsweise für f{x) = arcsin^ x (nämlich 
(x^ — 1)/"' + 3x/" + /' = 0), erzeugen. Doch damit nicht genug: Für die Koeffizienten 

der formalen Potenzreihe von aresin x = J2 CLnX^ folgt dann automatisch die holonome 

n=0 

Rekursionsgleichung n(l + n)(2 + n)a„+2 = n^^n, welche (glücklicherweise) nur die zwei Terme 
a„+2 und a„ enthält, sich daher lösen läßt und zu der Darstellung 



oo 4« „ 2 

• 2 ^ "'• 2n+2 

arcsm x = > — —x 



n=0 



+ n) (1 + 2n)! 



führt (vgl. [II, m, 0-0). 

Oder notieren wir uns lediglich die holonomen Rekursionen rationaler Funktionen sowie die 
der Funktionen 

(mn + b)\ , — (m G Z) sowie (1) 

(mn + by. 

(d. h. wir geben sie unserem favorisierten Computeralgebrasystem in geeigneter Form bekannt), 
so lassen sich nun für alle möglichen durch Addition und Multiplikation erzeugten Funktionen 
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holonome Rekursionen herleiten, beispielsweise die beiden Rekursionen 

{n-k + ifF^n + 1, fc) - (1 + nfF{n, k) = (2) 

und 

{k + lfF{n, k + 1)- (n- kfF{n, k) = (3) 

für F(n, k) = (j^ . Diese können zwar auch leicht direkt aus der Darstellung von F(n, k) 
abgelesen werden, aber das vorgestellte Verfahren kann problemlos auf die kompliziertere 
Funktion -F(n, k) = angewandt werden und führt dann auf die Rekursionen 

= nF{n + 2, /c) - (1 + 3n + n'^)F{n + 1, A;) + (1 + n^Fin, k) 

und 

Q = k{2 + kfF{n, k + 2)-{l + k){l + 3A; + k'^){?, + 3A; + k'^)F{n, k + l) + k{l + kfF{n, k) . 

Eine wichtige Fragestellung der Kombinatorik ist, zu einer gegebenen Funktion F{n,k), 
ausgedrückt als Produkt von Termen der Form (|l]), die unendliche Summe 

s{n) = ^ F{n, k) 

k 

zu berechnen, wobei über alle A; G 2 zu summieren ist. (In der Praxis sind jedoch häufig 
nur endlich viele Terme von Null verschieden.) Ist nun F{n,k) ein hypergeometrischer 
Term, d. h., sind sowohl F{n + l,k)/F{n,k) als auch F{n,k + 1)/F{n,k) rational bzgl. 

beider Variablen n und k - dies trifft beispielsweise wegen für F{n, k) = (j^ zu - 

so findet der (schnelle) Zeilberger- Algorithmus (@], s. auch und p3|) ^ine holonome 
Darstellung, d. h. eine holonome Rekursion, für s{n). Dieser Algorithmus baut auf dem von 
Gosper [14| gefundenen Entscheidungsalgorithmus für die unbestimmte Summation auf. In 

unserem Beispielfall findet der Zeilberger- Algorithmus für s{n) =!](?) = J2 w) die 

holonome Rekursion 

(1 + n) s{n + 1) = 2(1 + 2n) s{n) , 
welche (glücklicherweise) wieder nur zwei Terme hat. Daher bekommen wir die Darstellung 



s{n) = J2 

k=0 



Im allgemeinen wird die resultierende Rekursion bzw. Differentialgleichung natürlich mehr als 
zwei Terme enthalten. Aber auch dann enthält diese zum einen eine interessante Strukturinformation 
(z. B. über die Orthogonalität eines Polynomsystems) und kann zudem eine für numerische 
Zwecke nützliche Vorschrift darstellen (vgl. ||TÖ[|-|[TT[|). 

Die gefundene Strukturinformation kann insbesondere zur Identifikation transzendenter 
Funktionen herangezogen werden. Um z. B. die Identität 



fc=0 \'^/ k=0 



n 
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(vgl. p8|) zu überprüfen - diese ist nicht trivial, da die beiden Summanden nicht dieselben 
Rekursionen bzgl. k und n erfüllen! - brauchen wir nur zu zeigen, daß beide Summen derselben 
Rekursion 

{n + 2ys{n + 2) - (16 + 21n + 7n^)s{n + 1) - {n + l)^s(n) = 

genügen - dies macht der Zeilberger- Algorithmus - sowie dieselben Anfangswerte s(0) = 1 
und s(l) = 2 haben (dies ist trivial). 

Als weiteres Beispiel betrachten wir die Funktion (a, /?, 7 G INq) 

V( ß ) = r_ir+/5+^ T{a+ß+j-d/2)Tid/2-^)Tia+j-d/2)Tiß+j-d/2) 
^ '^'^^ ^ ^ T{a)T{ß)T{d/2)T{a + ß + 2^-d){m^)^+ß+^'^ 

/ a + ß + -f - d , a + 7-d/2 
a + ß + 2-f-d 

{2F1 stellt hierbei die Gaußsche hypergeometrische Reihe dar, s. z. B. Kapitel 15), welche 
bei der Berechnung von Feynman-Diagrammen eine Rolle spielt |T2|, für die wir die Rekursion 

= {2a-d+2-f) {2a + 2ß-d+2-f) {2 + 2a + 2ß-d+2-f) V{a,ß,-f) 
-2a {2 + 2a + 2ß-d+2-i)m^ 

■ {-2a-2ß+2d- A-f + 2z+Aaz + 2ßz-3dz+A-fz) V{l + a, ß, 7) 
+8 a {1 + a) {1 + a + ß - d + {z - 1) zV{2 + a, ß, 7) 

sowie analoge Rekursionen bzgl. der Variablen ß und 7 erhalten. Diese können dann zur 
numerischen Berechnung herangezogen werden. 

Zeilberger betrachtete in |^ die allgemeinere Situation von Funktionen F mehrerer diskreter 
und stetiger Variabler. Sind es d Variablen und hat man d (im wesentlichen unabhängige) 
möglicherweise gemischte homogene partielle Differential-Differenzengleichungen mit Polynomkoeffizienten 
(bzgl. aller Variablen) für F, nennen wir F ein holonomes System (vgl. [§-[0|). Dann legen 
diese Gleichungen F zusammen mit geeigneten Anfangswerten bereits eindeutig fest. 
Insbesondere gilt dies also, wenn das gegebene System holonomer Gleichungen separiert ist, 
d. h., wenn in jeder der Gleichungen nur Ableitungen bzgl. einer der stetigen Variablen bzw. 
nur Shifts bzgl. einer der diskreten Variablen vorkommen. Beispielsweise bilden die Legendre- 
Polynome F{n,x) = Pn{x) (|I[], Kapitel 22) auf Grund ihrer Differentialgleichung 
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[X 



l)F"(n, x) + 2xF'{n, x) - n{l + n)F{n, x) = (4) 



sowie ihrer Rekursionsgleichung 

{n + 2)F{n + 2, x) - (3 + 2n)xF{n + 1, x) + (n + l)F(n, x) = (5) 

zusammen mit den Anfangsbedingungen 

F(0,0) = 1, F(1,0) = 0, F'(0,0)=0, F'(1,0) = 1 

ein holonomes System. Gleichungen (^) und (^ stellen also eine hinreichende algebraische, ja 
sogar polynomiale, Struktur zur Beschreibung von -P„(x) dar. 
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Faßt man die auftretenden (partiellen) Differentiationen und Indexverschiebungen als Operatoren 
und die Differenzen-Differentialgleichungen als Operatorengleichungen auf, so stellen diese 
ein polynomiales Gleichungssystem in einem nichtkommutativen Polynomring dar. Ist 
nämlich x eine stetige Variable und der zugehörige Ableitungsoperator, so ist wegen der 
Produktregel D.j.{xf) — xD^f = /, und folglich hat man den Kommutator D^x — xD^ = 1. 
Ist andererseits k eine diskrete Variable und K der zugehörige Indexverschiebungsoperator 

= ük+i ( Aufwärt sshift), so gilt K{kak) — kKük = {k + l)ak+i — kau+i = «fc+i = 
und folglich gilt die Kommutatorregel Kk — kK = K. Entsprechendes gilt für die restlichen 
Variablen, während alle anderen Kommutatoren verschwinden. 

Das Umformen eines durch gemischte Differenzen-Differentialgleichungen gegebenen holonomen 
Systems stellt sich in dem betrachteten nichtkommutativen Polynomring als ein polynomiales 
Eliminationsproblem dar, welches mit nichtkommutativen Gröbnerbasen-Methoden gelöst 
werden kann (H, [|l|], pf, [||], ü, [ü-il). Als Beispiel diene F{n,k) = Hierfür 



gilt die Pascalsche Dreiecksbeziehung F{n + 1, k + 1) = F{n, k) + F{n, k + 1) sowie die reine 
Rekursion (n + 1 — k)F{n + 1, k) — {n + l)F{n, k) = bzgl. n, welche in Operatornotation 
{KN-l-K)F{n, k) = sowie ({n+l-k)N- (n + l))F(n, A;) = lauten, wobei NF{n, k) = 
F{n+1, k) den Verschiebungsoperator bzgl. n bezeichne. Somit erhält man das Polynomsystem 

KN - 1 - = sowie (n + 1 - k)N - (n + 1) = . 

Die Gröbnerbasis des erzeugten Linksideals bzgl. der Termordnung (/c, n, N) (lexikographisch) 



ergibt sich mit ||2^ zu 

[{k + l)K + k - n,{n + 1 - k)N - {n + 1), KN -1-K] 
d. h. also, daß auf diesem Wege automatisch die reine Rekursion 



{k + 1) F(n, k+l) + {k-n) F{n, fc) = 



bzgl. k erzeugt wurde. 

Als weiteres Beispiel betrachte ich die Legendre- Polynome, für die die Beziehungen (@)-(^ 
gelten. Hier haben wir also das Polynomsystem [D := D^) 

(x^ - 1)0"^ + 2xD - n(l + n) = sowie {n + 2)N^ - (3 + 2n)xN + (n + 1) = . 
Die Gröbnerbasis des erzeugten Linksideals bzgl. der Termordnung (D, A^, n, x) ergibt sich zu 

[{x^ - 1)D^ + 2xD - n(l + n), 
(1 + n)ND - (1 + n)xD - (1 + nf, 

{x^ -1)ND - {l + n)xN + (l+n), (6) 
{l + n){x^-l)D-{l + nfN + x{l + nf (7) 
(n + 2)N^ - (3 + 2n)xN + (n + 1)} , 
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wobei ich der besseren Lesbarkeit halber die Operatoren D und N wieder rechts positioniert 
habe. Hier wurden also D-Potenzen weitestgehend eliminiert, und Gleichungen entsprechen 
den Beziehungen 

ix' - l)PUiix) = il+n) {xPn+i{x) - P„(x)) (8) 

(x^ - l)P;(x) = {l + n) (P„,+i(x) - xP„(x)) (9) 

zwischen den Legendre-Polynomen und ihren ersten Ableitungen. Man sieht also, daß auf 
diesem Wege neue Beziehungen (zwischen den Binomialkoeffizienten bzw. zwischen den Ableitungen 
der Legendre-Polynome) hergeleitet wurden. 

Analog lassen sich mit dieser Methode Rekursionen für holonome Summen herleiten. 
Betrachten wir beispielsweise 

n n / \ 

s{n) = Y.nn.k) = Y.[APn{x), 

k=0 k=0 V^/ 

dann findet man mit dem Produktalgorithmus zunächst die holonomen Rekursionen 

{n-k + l)F{n + 1, Ä;) - (1 + n)F{n, k) = 

sowie 

(2 + kYF{n, k + 2) ~ {3 + 2k){n - k - l)xF{n, k + 1) + {n - k){n - k - l)F{n, k) = 
für den Summanden F{n,k). In der Gröbnerbasis des von den zugehörigen Polynomen 

(n- A; + 1)A^- (1 + n) 

sowie 

(2 + kfK' - (3 + 2k){n -k- l)xK + (n - k){n - k - l) 
bzgl. der Termordnung (/c, n, N) erzeugten Linksideals liegt das fc-freie Polynom 

(2 + nyX^N' - K{2 + n){3 + 2n){K + x)N + (1 + n)(2 + n){l + K' + 2Kx) , 

welches einer /c-freien Rekursion für P(n, k) entspricht. Da bei der Summation über G 2 
die verschobenen Summen 

s{n) = J2 F{n, k)=J2 ^(^. ^ + 1) = E ^(^' ^ + 2) 

k k k 

alle dieselbe Summenfunktion s{n) liefern, liefert die Substitution K := 1 die gültige holonome 
Rekursion 

(2 + n)s{n + 2) - (3 + 2n)(l + x)s{n + 1) + 2(1 + n)(l + x)s{n) = 

für s{n). 

Die gezeigte Methode ist zur Generierung von Identitäten im Prinzip universell einsetzbar, 
benötigt jedoch den komplizierten Apparat des (nichtkommutativen) Buchberger- Algorithmus 
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und erbt die damit verbundenen Nachteile. Eine wesentliche Hürde stellt die Komplexität bei 
Problemen mit vielen Variablen dar. 

Interessiert man sich wie bei obigem Beispiel für die Erzeugung von Identitäten 

zwischen den Ableitungen F^^\n + fc, x) (j, G Nq) eines holonomen Systems F(n,x), 
muß dies aber nicht unbedingt sein. Es geht in vielen Fällen bereits mit linearer Algebra! 
Dazu muß man aber mehr Information hineinstecken. Die geeignete über die holonomen 
Beziehungen hinausgehende Information besteht in einer Beziehung der Form 



m 

F'{n, x) = ^ rk{n, x)F{n + A;, x) 

k=0 



mit rationalen Funktionen bzgl. n und x, einer Ableitungsregel für F{n,x) also, sofern 
erhältlich. Es zeigt sich, daß in der Praxis holonome Systeme (wie beispielsweise Systeme 
orthogonaler Polynome etc., s. z. B. 22.8, und JlÖl) so strukturstark sind, daß eine 
Ableitungsregel verfügbar ist. Man kann zeigen, daß Summe und Produkt solcher Systeme in 
der Regel auch wieder holonome Systeme mit Ableitungsregel darstellen , und Abhängigkeiten 
zwischen den Ableitungen F'^^\n + k^x) {j,k G Nq) können mit reiner linearer Algebra 
gefunden werden. 

Auf diese Weise wurde z. B. die Beziehung 

p{a,ß)(^^ ^ 2{a + n){ß + n) p("'/^)Vr^ 

" ^ ' {a + ß + n){a + ß + 2n){a + ß + 2n+l) "^^ 



{a + ß + 2n){a + ß + 2n + 2) 



x] 



I 2{a + ß + n+l) pi-ßYi^) 

{a + ß + 2n + l){a + ß + 2n + 2) 

für die Jacobi- Polynome Pj^^'^^x) Kapitel 22) automatisch erzeugt. Hierbei war die 

Zielsetzung, daß die auftretenden Koeffizientenfunktionen von -P^^'f ■* {x) nicht von x abhängen 
sollen. Dies ist für Fragestellungen aus der Spektralapproximation (s. 0, § 2.3.2) von Bedeutung 
Zuletzt will ich darauf verweisen, daß die vorliegende Beschreibung natürlich keinen Anspruch 
auf Vollständigkeit erheben kann. Ich konnte weder auf den Gosper- Algorithmus noch 



auf die Wilf-Zeilberger-Theorie der WZ-Paare und rationalen Zertifikate eingehen [p3| -p7 



Auch Petkovseks Algorithmus p^, der alle hypergeometrischen Termlösungen holonomer 
Rekursionsgleichungen berechnet, konnte keine Berücksichtigung finden. 
Die in diesem Artikel durchgeführten Rechnungen wurden mit Mathematica und Reduce 
durchgeführt. Mein Dank gilt Prof. Peter Deuflhard, der mich ermutigt hat, mich mit dem 
vorliegenden Thema zu beschäftigen, Herbert Melenk, mit dem ich wichtige Gespräche über 
nichtkommutative Gröbnerbasen führen konnte, sowie Jochen Fröhlich, der mich auf die 
Spektralapproximation aufmerksam machte. 
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